Quantization of Gauge Theory for Gauge Dependent Operators 
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Based on a canonically derived path integral formalism, we demonstrate that the perturbative 
calculation of the matrix element for gauge dependent operators has crucial difference from that for 
gauge invariant ones. For a gauge dependent operator 0(0) what appears in the Feynman diagrams 
is not 0{(j>) itself, but the gauge-transformed one OC^fj})^ where ui characterizes the specific gauge 
transformation which brings any field variable into the particular gauge which we have adopted to 
quantize the gauge theory using the canonical method. The study of the matrix element of gauge 
dependent operators also reveals that the formal path integral formalism for gauge theory is not 
always reliable. 
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Ever since the emergence of nucleon spin problem, the 
matrix element of gauge dependent operators has aroused 
great interest ||^ . This is because the conventional under- 
standing of gluon spin, gluon and quark orbital angular 
momentum all corresponds to gauge dependent opera- 
tors. When studying these operators, people have natu- 
rally adopted the same calculation scheme as for gauge 
invariant ones, however some quite contradictory results 
were obtained In this paper, by quantizing the gauge 
theory from the very beginning, we demonstrate that the 
perturbative calculation rules for gauge dependent oper- 
ators have nevertheless critical difference from those for 
gauge invariant ones. By doing so we also reveal that the 
highly formal path-integral formalism for gauge theory is 
not always reliable. 

The theory we start with is the classical SU(N) or U(l) 
gauge theory. We adopt the well-defined canonical quan- 
tization approach to quantize it. Due to gauge freedom, 
we must specialize a gauge at first. Here we choose the 
most convenient axial gauge 



n • ^ = 0, 



(1) 



where A is the gauge potential, n is an arbitrary but fixed 
constant four-vector. Now we can apply the standard 
canonical quantization procedure and obtain the com- 
mutation relations and the quantum Hamiltonian. But 
instead of going further to use the operator method to 
derive the Feynman rules, we use the so obtained quan- 
tum Hamiltonian to derive the path-integral expression 
for the vacuum expectation value of an arbitrary operator 
0{(j)), where (p denotes collectively the gauge and mat- 
ter fields. After the standard procedure which can be 
found in textbooks we can write up to an irrelevant 
normalization factor 



The notations are standard. This expression is rigor- 
ous and applies to both Abelian and non-Abclian gauge 
theory Q|. Moreover, we have no requirement for the 
operator O. It can be either gauge dependent or gauge 
invariant. However, Eq. (^) is non-perturbative and not 
yet perturbatively applicable. To derive the rules for per- 
turbative calculation, we use the Faddeev-Popov trick by 
multiplying Eq. (H) with the identities 



1 = Af{A, x) J [Diu]6iFi-A) ~ x), (3a) 
1 = / [Dx]G[x], (3b) 



where F and G are two arbitrary functions (or in general 
functionals), [Duj] is the invariant measure of the gauge 
group, '^A is the result of A after a gauge transformation 
characterized by w. F is usually called the gauge fix- 
ing function, for example we can take our original axial 
gauge form F{A) = n-A. Accidentally, it is apparent that 
by this method we can transit to an arbitrary gauge by 
choosing arbitrary F and G (this point will be explained 
further later). 

Multiplying (||) with (^) and (3b), we get 



= J m 



e'^C>((/))<5(7i • A). 



(2) 



Z{0)= J [D(j)][Duj][Dxy^OS{n ■ A) 

xAf{A,x)S{F{'^A)-x)G[x]. (4) 
Now make the transformation 

0^-^10, (5) 

note that [Dfjj], S, and Ap{A, x) are gauge invariant, and 
integrate over x, we get 
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Z{0)^ J[Dcj,y^AF{A,F{A)G{FiA)) 

X / [Duj]{'^'' 0)5{n A) (6) 



The second line in Eq. gives up to an irrelevant 
constant factor [|j 

""O, (7) 
where ujq is the specific function of A which satisfies 



n-'^° A = 0. 
So we obtain the aimed expression 



(8) 



Z{0)= [D(j)]e''^{'^'>0)AFiA,F{A))G{F{A)). (9) 



By choosing a Gaussian type of G, and introduce the 
ghost field to treat the Faddeev-Popov determinant A^?, 
one can straightforwardly derive Feynman rules from Eq. 

®- . 

We note very importantly that what appears in this 

final expression is '^"O instead of O itself. This will give 
critical difference for gauge invariant and non-invariant 
operators. 

First we look at the gauge invariant operator: '^"O 
equals O. Now Eq. (||) becomes 

Z{0)^ [ D(t)e'^OAFiA,F{A))GiF{A)). (10) 



So we see that starting from a canonical quantization in 
the specific axial gauge, we finally arrive at an expression 
independent of the original gauge condition! For instance 
the gauge A2 = or ^3 = lead to the same Eq. (10). 



This states that the Green's function of gauge invariant 
operators are gauge independent. 

For gauge dependent operators, however, '^°0 does not 
equal O. So the perturbative calculation of the Green's 
function for gauge dependent operators is not a straight- 
forward generalization of that for gauge invariant ones. 
Instead, we encounter an additional complication of al- 
ways having to apply a gauge transformation to the op- 
erator before inserting it into the Feynman diagrams. 

It should be remarked that in Eq. (||) only the op- 
erator O suffers such a gauge transformation, the other 
field variables such as those in the action remain the ordi- 
nary form. Accordingly, in the perturbative calculations, 
only the external vertex corresponding to O should be 
the gauge transformed one, while the Feynman rules for 
the internal vertices and propagators generated from the 
expansion of the effective action are still the ordinary 
ones. 

We note that in Eq. (^) ujo is universally given by Eq. 
(^. However, the form of F and G can be arbitrary, i.e. 
the Feynman rules can be arbitrary but they give the 
same final results. One may wonder why Eq. (H) takes 



such an apparently peculiar form, especially why must 
we apply a transformation to O and why loq still takes 
the form in Eq. (^) even after we have "transitted to 
another gauge" by choosing an F(A) different from n ■ A. 
We explain below that this actually has deeper physical 
reasons behind. 

Due to the extra degree of gauge freedom, we have to 
fix a gauge to do quantization. And such a gauge choice 
will naturally leave "traces" on the considered operator 
O and also on the subsequent formalism, such as the 
to-be-derived quantum Hamiltonian and Feynman rules. 
Indeed, different gauges lead to different form of quantum 
Hamiltonian. Or to say, the form of quantum Hamilto- 
nian is gauge dependent. Furthermore, if we stick to the 
operator method to calculate the Green's functions, the 
original gauge choice will go throughout to manifest in 
the Feynman rules. However the gauge invariance prin- 
ciple commands that we must in the end show all the 
physical observables to be independent of the original 
gauge we chose in quantization. Within the canonical 
operator method, such a demonstration is essentially a 
direct proof of the equivalence of Feynman rules derived 
in different gauges, and will necessarily involve us in a de- 
tailed analysis of the Feynman diagrams. This was done 
for QED by Feynman H and for Yang-Mills theory by 
Cheng and Tsai @] . 

In contrast, instead of demonstrating gauge invari- 
ance at the level of Feynman diagrams, the advantage 
of path-integral formalism is that we can already get rid 
of the original gauge condition at an intermediate stage 
before deriving the Feynman rules. (To be equivalent to 
the canonical formalism, our path-integral expression is 
derived by the operator method using the canonically- 
derived quantum Hamiltonian.) The crucial step is that 
we can obtain Eq. (|^), in which what appears is the 
original gauge invariant classical action S. (This can 
be done at least in the axial gauge, in which the quan- 
tum Hamiltonian takes a relatively simpler form, than, 
say, in the Coulomb gauge.) Now, the subsequent pro- 
cedure in Eqs. (0)-(||) is essentially converting n ■ A into 
n ■ A' — n ■'^ A — {F{A) — x) by coordinate transfor- 
mation, then integrating over x with a weight-functional 
G[x]- The Af{A,x) is nothing but the Jacobian of this 
coordinate transformation, which due to the compact- 
ness of the gauge group can be chosen gauge invariant. 
This step is mathematically trivial, and by it alone we 
are still in the original axial gauge, but now expressed 
a.s n ■ A' = 0. However, since [D'^cj)], ApC^Ajx), and 
critically the action S{"(f>) are gauge invariant, we can 
changed their '^0 back to (j). Now we are in another ar- 
bitrarily chosen gauge {F{A) — x = 0), and can commit 
an arbitrary weight G[x] to each x- Therefore, the free- 
dom to choose arbitrary F and G is a consequence of 
the gauge invariance of S in Eq. (||). Namely, despite 
the appearance of an explicit S{n ■ A), the original gauge 
condition n ■ A has essentially disappeared out in Eq. (0) 
(apart from 0{(j))), because we can convert n ■ A into an 
arbitrary form by changing A to '^A but leave the A else- 
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where unaltered. Therefore later on when we convert G 
and Ai? into an effective action and expand it to derive 
Feynman rules, we already know for sure that these Feyn- 
man rules are equivalent to each other, without having 
to study the details of Feynman diagrams. This is the 
convenient path-integral method of demonstration that 
gauge theory still reserves gauge invariance after quanti- 
zation. Namely, the "trace" of the original gauge choice 
has disappeared out in the final Feynman rules. 

However, the "trace" left on the considered operator 
O by the original gauge condition will not disappear, if 
O is not gauge invariant. And such a "trace" should be 
"visible" no matter whatever tricks we play. Indeed, in 
the original expression Eq. the gauge condition is 
enforced by the delta function. And after we adopted 
the Faddeev-Popov trick to transit to an arbitrary gauge 
fixing form G{F{A)), this G{F(A)) no longer enforce the 
original gauge condition on O M , however we get a gauge 
transformation on O characterized by loq. In fact, by the 
definition in Eq. (||), the role of ujq on O is nothing but 
bringing it to the initial axial gauge in which we start the 
canonical quantization. 

In one word, once we start with the expression Eq. 
(^, which is obtained by canonical quantization in the 
axial gauge, the subsequent calculation will exhibit such 
a gauge condition anyhow, unless the considered operator 
O is gauge invariant. And no matter what kind of gauge 
fixing functions F and G we choose to derive the Feyn- 
man rules, the final matrix element we obtain should be 
the same (as long as the Faddeev-Popov trick is justified) 
and should be regarded as the result in the original axial 
gauge; the freedom of choosing Feynman rules for calcu- 
lation is simply because the quantized gauge theory still 
reserves gauge invariance, as was demonstrated by the 
path-integral formalism. If we want to study what the 
matrix element of a gauge dependent operator would be 
in another gauge, we have to begin the initial canonical 
quantization in that gauge (which might not necessarily 
lead to a simple form like Eq. Merely adopting the 

Faddeev-Popov trick to shift the form of F and G only 
changes the gauge for Feynman rules (which are equiva- 
lent in different gauges), but not the gauge for the whole 
matrix element. 

For example, the conventional understanding of quark 
or electron orbital angular momentum corresponds to the 
gauge dependent operator 

CpX%l>^ X -i^ V. (11) 

Whether Lq gives gauge invariant matrix clement in a nu- 
cleon helicity eigenstate will determine whether quark's 
orbital contribution to nucleon spin is meaningful. Ac- 
cording to our above discussion, we cannot calculate i^'s 
matrix element as we do for the gauge invariant opera- 
tors such as quark spin. Actually we don't really know at 
all the expression of LqS matrix element in an arbitrary 



gauge except for the gauges in which we know how to do 
canonical quantization. 

In the usual discussion of gauge theory with path- 
integral approach, people sometimes simply ignore the 
gauge freedom and write down a formal expression as if 
all field variables were independent: 

Z{0) = j [Dcj.y'^Oicf). (12) 

Then by following the same steps as from Eqs. (|^) to (|^), 
we can get 

Z{0)= J[Dcj)]e'^Ap{A,F{A)G{F{A)) 

X [[Du^r-'O). (13) 



Different from Eq. (^, we don't have delta function in 
the second line to pick up a specific w as a function of A. 
The uj here is purely an integration variable independent 
of the field variable. If O is gauge invariant and hence 
O — O, the infinite integration J[Duj] factorizes out 
as an irrelevant normalization factor, so we get the same 
perturbatively applicable expression as what we obtained 
by starting with strict canonical quantization. 



Z{0)= [D(l)]e'''OAF{A,F{A)G{F{A)). 



(14) 



Since Eq. ( p^ can be derived from the same Eq. ( [T^ ) 
for arbitrary F and G, this formally states that Z{0) is 
independent of F and G up to an irrelevant constant fac- 
tor, or to say, the Feynman rules in different gauges are 
equivalent. This is again the same conclusion as wc dis- 
covered by starting with well-defined canonical method. 

However, we remark that such a formal procedure of 
starting with the ill-defined infinite expression Eq. ( |T^ ) 
is not always reliable. Such unreliability simply often 
does not manifest for gauge invariant operator, but can 
be seen by studying gauge dependent operators. Gener- 
ally speaking, when O is gauge dependent we are unable 
to arrive at the perturbatively applicable expression Eq. 
( p^ ) from Eq. (^). But there are some exceptions. For 
instance we take the product of electron field at different 
points: 



'O wiU be 



(15) 



(16) 



Therefore when x y, i/j'' {x)ip{y) is gauge dependent. 
However, we see in Eq. (^ that the gauge transformed 
phase factor e*"^^^""^'^*'^ can be absorbed into the irrel- 
evant integration over lu. Therefore we can obtain the 
same expression for Z{tfj'' {x)tp{y)) as for gauge invari- 
ant operators, and would conclude that 2'(-0'''(x)V'(y)) 
can be calculated with the usual Feynman rules in an 
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arbitrary gauge and is gauge independent! This however 
contradicts the results obtained by starting with the well- 
defined canonical quantization. And in fact, we known 
that the Fermion two point function is gauge dependent. 
A toy model demonstration of how such error arises can 
be found in Rcf. §. 

By some more delicate tricks, the authors in Ref. 
10 adopted the above formal path-integral formalism to 
demonstrate that if we use the same calculation rules as 
for gauge invariant operators, the gauge dependent quark 
orbital angular momentum operator will give gauge inde- 
pendent matrix element for a nucleon helicity state. This 
however is refuted by explicit 1-loop calculations . 

Finally, we mention that Cheng and Tsai has 
ever pointed out that the Feynman rules derived via the 
Fadeev-Popov approach often exhibit singularities (for 
example in the Coulomb and axial gauges). Such singu- 
larities must be dealt with care, otherwise at the two-loop 
level of non-Abelian gauge theory the Feynman rules in 
the Coulomb, covariant, and axial gauges derived via the 
Fadeev-Popov approach will be inconsistent with each 
other 0. 

In summary, we demonstrated that the calculation of 
the matrix element of gauge dependent operators is not a 
straightforward generalization of that for gauge invariant 
ones, and must be carefully derived from the very begin- 
ning of gauge theory quantization. We must distinguish 
the effects of initial gauge choice (in quantization) on the 
Feynman rules and on the studied operator or matrix 
element. Feynman rules in different gauge are demon- 
strated to be equivalent but the gauge condition on the 
studied operator might bring non-trivial difference. The 
formal path-integral formalism of starting with an ill- 
defined infinite expression, which often gives correct re- 
sults for gauge invariant operators, is however not always 
reliable, as can be seen by studying the matrix element 
of gauge dependent operators. 

Note added: We call the interested readers' attention 
to a most recent paper by Hoodbhoy and Ji . We got 
quite some inspiration from [ p^ in developing this paper, 
however our opinions towards gauge theory quantization 
are rather different from that in |13|. 
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